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Abstract. We report the measurement of charged hadron production at different 
pseudo-rapidity values in deuteron+gold as well as proton-proton collisions at ^snn = 
200GeV at RHIC. The nuclear modification factors RdAu and Rep are used to 
investigate new behaviors in the deuteron-|-gold system as function of rapidity and 
the centrality of the collisions respectively. 



1. Introduction 

The work reported in this presentation was motivated by the possibility of observing the 
onset of gluon saturation in d-l-Au collisions at RHIC. Saturation of the parton densities 
as their fractional longitudinal momenta x tends toward zero is widely accepted on the 
basis of hard limits set by unitarity (Froissart bound) P, but the proper scale for that 
new regime, as well as its properties are the subject of intense study. Deep inelastic 
scattering of leptons on hadronic systems performed at HERA and FNAL have shown 
that for small-x values gluons are the dominant component of the hadron wave function 
and that their densities grow as powers of ^ I2j. Limits to that growth have been 
proposed in a phenomenological description of the cross section of the virtual photon 
exchanged in those reactions. A functional form of that cross section is proposed such 
that it grows to become a constant beyond a scale that depends on a; as ~ where 
A = 0.2 — 0.3 is the only free parameter with values extracted from fits to data [3 E]. 
A QCD based theory of dense partonic systems accessible with the new high energy 
colliders has been developed And the new regime it strives to describe is called the 
Color Glass Condensate (CGC). This theory offers a novel way of calculating nuclear 
properties at high energies based on classical field theory techniques. 

The formalism of saturation can also describe p-|-A systems and the characteristic 

scale in transverse momentum is directly proportional to the density of partons in 

1 

the ion: ~ ^ and recalling the relation between the rapidity of a parton that 
carries the fraction x of longitudinal momentum and the rapidity of the hadron: 
Vparton = Uhadron — ^^(^) 0^16 cau Write ~ A^e^^ . RHIC offcrs the advantage of 
combining high energy and high A ions (197 for gold) and BRAHMS, one of the four 
experiments currently collecting data at RHIC, can study particle production at high 
rapidity. 
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2. Experimental setup 

The data presented in this report were collected with both BRAHMS spectrometers, 
the mid-rapidity spectrometer (MRS) and the front section of the forward spectrometer 
(FFS), complemented with an event characterization system used to determine the 
geometry of the collisions. A detailed description of the BRAHMS experimental setup 
can be found in 0. However, the low multiplicity of charged particles in the proton or 
light ion reactions required an extension of the basic apparatus with a set of scintillator 
counters (called INEL detectors). These detectors [TU] cover pseudo-rapidities in the 
range: 3.1 <| r] \< 5.29, and define a minimum biased trigger. This trigger is estimated 
to select ~ 91% ±3% of the 2.4 barns d-l-Au inelastic cross section and 71% ±5% of the 
total inelastic proton-proton cross section of 41 mb. The INEL detector was also used 
to select events with collision vertex within ±15 cm of the nominal collision point with 
a resolution of 5 cm. 

The centrality of the collision was extracted from the multiplicity of the event 
measured within the angular region | r] |< 2.2 with a combination of silicon and 
scintillator counters [Zj. The results reported here were obtained without making use 
of the full complement of particle identification detectors. Ongoing analysis of the data 
includes particle identification and the use of the complete forward spectrometer FS. 

3. Spectra 

Experimental methods and analysis techniques used to extract the spectra presented 
here are described in |8j. Figure ^ shows the invariant yields obtained from p+p 
collisions (panel a) and d+Au collisions (panel b). For each system we studied particle 
production at 40 degrees with the MRS spectrometer and 12 and 4 degrees with the FFS 
spectrometer. Each distribution was obtained from several magnetic field settings and 
corrected for the spectrometer acceptance, tracking and trigger efficiency. No corrections 
were applied to the spectra for absorption or weak decays. Statistical errors are shown 
as vertical lines, and an overall systematic error of 15% is shown as shaded boxes on 
all points. The data points for the spectra extracted at ?7 = 3.2 as well of those from 
7] = 2.2 have been placed at the mean px calculated within each bin. This turns out to be 
equivalent in magnitude to the unfolding of the momentum resolution. The momentum 
resolution of the spectrometers at full magnetic field is 6p/p = 0.0077p for the MRS 
spectrometer and 6p/p = O.OOlSp for the FS with p given in units of GeV/c. The 
P+P spectra have also been corrected for trigger efficiency by 13 ± 5% to make them 
minimum biased with respect to the total inelastic cross section. The proton+proton 
data sample at forward angles was collected in its majority with the magnets set to 
favor the detection of negative particles. The spectra at r/ = 3.2 have been fitted with 
a power law function q and the integral of that function over is compared for 

consistency in table Q with UA5 results [9^ for the p+p system, as well as our own 
multiplicity measurement for the d+Au system |lUj . 
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Figure 1. Spectra for charged hadrons at different pseudo-rapidities. Panel a shows 
the spectra obtained from proton-proton collisions and panel b those from d-f-Au 
collisions. The top most distributions in both panels correspond to the invariant yields 
of measured at 40 degrees with the MRS spectrometer (scaled by 100 for clarity 

purposes), followed by the yields of negative hadrons measured at 12 (scaled up by 10) 
and 4 degrees respectively. The distributions obtained at 77 = 2.2 and 3.2 have variable 
bin size indicated by horizontal lines, the data points have been located at the mean 
value calculated within that particular bin. The distribution obtained at 77 = 1 has 
fixed bin width of 200 MeV/c and the data points are displayed at the middle of the 
bin. 



Table 1. Fits to power law shapes dX rj — 3.2. 



System 


dN 1 dN 

dv fit ' d.r] ^gas 


Pq 


n 


X^/NDF 


P + P 
d + Au 


1.05 ± 0.06 / 0.95 ±0.07 
2.23 ±0.09 / 2.1 ±0.6 


1.18 ±0.16 
1.52 ±0.1 


10.9 ±0.9 
12.3 ±0.5 


13. / 11 
102 / 11 



4. The nuclear modification factor 

The d-|-Au spectra are compared to normahzed p+p distributions measured at the same 
spectrometer setting. This comparison with p-|-p results is based in the assumption 
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that the production of moderately high transverse momentum particles scales with the 
number of binary collisions N^oii in the initial stages. The so called nuclear modification 
factor is defined as: 

1 NdAuiPT,'n) 



where Ncoii is estimated to be equal to 7.2 ± 0.3 for minimum bias collisions. 
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Figure 2. Nuclear modification factor for charged hadrons at pseudorapidities 
?7 = 0, 1.0, 2.2, 3.2. Statistical errors are shown with error bars. Systematic errors are 
shown with shaded boxes with widths set by the bin sizes. The shaded band around 
unity indicates the estimated error on the normalization to (Ncoii)- Dashed lines at 
Pt < I GeV/c show the normalized charged particle density ratio (jv^^ '^^N/d^(pp)'^ ■ 



Figure |21 shows the nuclear modification factor defined above for four rj values. 
The ratio at ?7 = and 1 was calculated by averaging on both positive and negative 
charges {h^ + h~)/2, and the one at rj = 2.2 and 3.2 has been calculated including 
only h~ particles. The average charge is the preferred quantity because it compensates 
the isospin mismatches between the systems being compared. Because we do not have 
enough data to extract the distribution of positive charged particles from p+p, the ratio 
at T] = 2.2 and 3.2 was calculated with spectra of negative charged particles. We can 
use the model PYTHIA to state that if we were to make the ratio with the average 
charge, the ratios would be smaller than those reported in this work. At mid-rapidity 
(?7 = 0), the nuclear modification factor exceeds 1 for transverse momenta with values 
greater than 2 GeV/c in similar way as the measurements performed by Cronin at lower 
energies JT] . Such enhancement has been related to multiple scattering at the partonic 
level [ig. 
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Table 2. Npart and Ncoii values extracted from HIJING calculations for d+Au 
collisions. 



Centrality 




Npartid) 




Central 0-20% 


12.5 


1.96 


13.6 ±0.3 


Semi — central 30 — 50% 


7.36 


1.79 


7.9 ±0.4 


Peripheral 60 - 80% 


3.16 


1.39 


3.3 ±0.4 



It is seen that a shift of one unit of rapidity is enough to make the Cronin type 
enhancement disappear, and further increases in t] decrease even further the scaled yield 
of charged particles produced in d+Au collisions. 

In all four panels, the statistical errors, shown as error bars (vertical lines), are 
dominant specially in our most forward measurements. Systematic errors are shown as 
shaded rectangles. An additional systematic error is introduced in the calculation of 
the number of collisions Ncoii that scales the d+Au yields to a nucleon-nucleon system. 
That error is shown clS 8b 15% band at RdAu = 1- 

We see a one to one correspondence between the RdAu values at low px and the 
ratio , as demonstrated in Fig. |21 where that ratio is shown as dashed 

{Ncoll) dN/dri(pp) '-' 

lines at Pt < 1- 



5. Centrality dependence 

The relation between particle yields and the centrality of the d+Au collision has 
also been explored using the fact that the multiplicity of charged particles in a 
reasonably wide t] range (| r] \< 2.2) correlates well to the impact parameter of 
the collision. Three samples of different centralities were defined according to the 
multiplicity of each event, and scaled histograms in transverse momentum were filled: 
NcentraiiPT) = Tr~ No-2o%(pt) foT cvcuts with multiplicities ranging from to 20%. 

coll 

semi— central (Pt) = ]v~" ^30-50% (Pt) for semi-central events with multiplicities ranging 
from 30 to 50%, and finally, Nperiph{PT) = ]v^^60-80%(pt) for peripheral events with 
multiplicities ranging from 60 to 80% with Ncou values listed in table El 

With these histograms, two ratios were constructed: i2??p*^'*' = ^^<^^trai{pT) ^^^^^ 

^semi— central _ Ng^^i_^^„tj.^i{pT) 

Because these ratios are constructed with events from the same data run, many 
corrections cancel out. The only correction that was applied to these ratios is related 
to trigger inefficiencies that become important in peripheral events. The dominant 
systematic error in these ratios stems from the determination of the average number of 
binary collisions in each centrality data sample. This error is shown as a shaded band 
around 1 in Fig. 01 

The four panels of figure El show the central (filled symbols) and semi- 

central /^^™-ce«t'-a-^ (open symbols) ratios for the for rj settings. Starting on the left 
panel corresponding to ?7 = 0, the central events yields are systematically higher than 
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Figure 3. Central (full points) and semi-central (open points) Rep ratios (see text for 
details) at pseudorapidities g = 0,1.0,2.2,3.2. Systematic errors (~ 5%) are smaller 
than the symbols. The ratios at all pseudorapidities are calculated for the average 
charge . 



those of the semi-central events. On the second panel corresponding to events measured 
at ?7 = 1, the scaled yields of both samples are almost equal. At ?7 = 2 the trend has 
already reversed, this time the semi-central events have higher yields than the central 
events, and finally, at = 3, the yields of central events are ~ 60% lower than the 
semi-central events for all values of transverse momenta. 

6. Discussion 

The first results from the 2003 RHIC run showed that Cronin like enhancements are 
still prominent around mid-rapidity at RHIC energies and that jet suppression detected 
in Au-Au systems [12] is not present in d+Au [HI jH]. It was then concluded that jet 
suppression must happen in a medium formed after the interaction of the Au ions, and 
that medium is partonic and highly opaque. 

In recent months, several groups reported theoretical discussions as well as 
numerical calculations within the context of the CGC with predictions about trends in 
the physics of deuteron-gold collisions at RHIC Several authors showed that the color 
glass condensate exhibits a Cronin like enhancement at mid-rapidity [TK t lTH | IT7 | ni^ while 
others argued for a gradual disappearance of such an enhancement as the energy of the 
collision or the rapidity of the measured particles is increased ^H] or offered numerical 
calculations based on similar theoretical arguments [201- All these theoretical works at 
higher rapidities introduce a modification of the gluon densities related to gluon emission 
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probabilities that grow as the rapidity gap between the probe and the target increases, 
within a formahsm called 'quantum evolution'. All these works predicted behaviors of 
the nuclear modification factor and its centrality dependence that are consistent with 
the measurements being reported in this paper. 

However, other groups based their predictions and the description of the data with 
a two component model that includes a parametrization of perturbative QCD and string 
breaking as a mechanism to account for coherent soft particle production, the so called 
HIJING model jSIj. Multiple scattering at the parton level has been added to this model 
with different degrees of sophistication and it describes well the d+Au results at mid- 
rapidity as well as the multiplicity densities that are dominated by soft processes. These 
models produce Cronin type enhancements at higher rapidities, and the magnitude of 
that enhancement grows with centrality. j22], j^Hj. Such behavior is not present in the 
data presented here. 

More recently Miklos Gyulassy pointed the higher yield of positive charged particles 
dX 7] = 3.2 shown in Fig |3] as 'clear evidence of valence quark fragmentation' [25 and 
the fact that the HIJING model provides a good representation of the data with a 
mechanism of string breaking introduced to explain low energy p-A measurements |25j . 
The authors of that paper offer an explanation for the triangular shape of the ratio 
'^^N/frii^p) ^^^^ present in these measurements. 



CM 

■D 



10 
1 

■ 10 ' 



'^-'-lO' 



T3 



^1 ^ 10 



10 



10 



10 



* negatives 

★ Positives 



3 4 
Pt [GeV/c] 



Figure 4. Spectra of positive and negative charged particles aX r/ — 3.2 

In summary, BRAHMS has studied the production of charged particles and 
compared their yields to those produced in proton-proton collisions scaled by the number 
of binary collisions and finds a clear evolution as the t] of the particles changes from 
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mid-rapidity where a Cronin type enhancement is seen above transverse momentum of 
2 GeV/c, to a gradual suppression as the values of rj change to 1, 2.2 and 3.2. We have 
also found that the suppression increases with the centrality of the collision. 
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